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In  the  past  several  years,  a numerical  method  has  been  developed  which  is 
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are  further  reduced  to  a simplest  possible  case  and  yet  still  contain  the 
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statement  has  been  established  which  is  both  the  necessary  and  sufficient  con- 
dition for  the  coupled  differential  equations  and  some  general,  but  physical 
meaningful  boundary  conditions.  The  finite  element  matrix  equations  are  then 
derived  from  this  variational  statement  Illustrating  the  way  that  coupling  terms 
could  be  handled  in  general. 
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motion  of  flapping  and  root  torsion  without  any  aerodynamic  force,  if  the  tor- 
sional spring  (or  the  pitch  control  link)  is  not  sufficiently  stiff.  This  in- 
stability does  not  appear  to  have  been  reported  previously. 
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1.  INTRODUCTION.  An  analytical  investigation  on  vibrations  and 
dynamic  stability  of  helicopter  rotor  blades  usually  consists  of  two  phases 
(1)  the  derivation  of  the  governing  differential  equations  to  include  para- 
meters and  variables  considered  physically  important,  and  (2)  the  formulation 
of  solutions  for  the  equations  derived  and  data  interpretations.  This  paper 
deals  with  the  second  phase  of  such  an  investigation. 

Due  to  the  slenderness  of  a helicopter  rotor  blades,  its  aerodynamic 
cross-section  and  the  requirements  on  the  craft's  maneuverability,  there 
are  a large  number  of  interacting  parameters  and  the  resulting  differential 
equations  are,  as  a rule,  nonlinear,  coupled  in  terms  of  field  variables. 

In  addition,  the  aerodynamic  forces,  coriolis  forces  due  to  rotation  of  the 
blade  are  nonconservative  in  nature  and  the  effects  due  to  structural  damp- 
ing and  various  boundary  conditions  must  be  evaluated. 

Considerable  attention  has  been  given  recently  to  the  derivation  of  a 
consistent  set  of  nonlinear  differential  equations  together  with  aero- 
dynamic forces.  This  fact  is  amply  demonstrated  by  the  work  of  Friedman 
and  Tong  [1]  and  that  of  Hodges  and  Ormiston  [2].  In  both  references  [1] 
and  [2],  brief  reviews  on  earlier  work  on  helicopter  blade  equations  can 
be  found.  As  for  obtaining  solutions  to  these  equations,  there  does  not 
appear  to  exist  a general  and  efficient  method  to  deal  with  the  difficulties 
associated  with  nonlinearities,  nonconservative  forces,  coupling  terms, 
various  boundary  conditions,  damping  effects  and  periodic  excitations.  For 
example.  Miller  and  Ellis  [3],  Ham  [4],  and  Friedman  and  Tong  [1]  have 
obtained  approximate  solution  by  including  in  their  solution  formulation 
only  the  lowest  modes  of  vibrations,  Hodges  and  Ormiston  [2,5]  employed  the 
Galerkin  technique  in  their  numerical  examples.  One  of  the  disadvantages 
of  this  approach  is  its  inability  to  handle  general  boundary  conditions. 

Using  the  combined  concept  of  adjoint  variable  and  Lagrange  multipliers, 
variational  statements  can  be  established  for  a wide  range  of  linear  prob- 
lems with  nonconservative  forces  and  very  versatile  boundary  conditions  [6]. 
Thus  a generalized  Rayleigh-Ritz  approximation  scheme  can  be  established 
for  the  obtaining  of  solutions  of  these  otherwise  difficult-to-solve  prob- 
lems. In  conjunction  with  finite  element  discretization,  this  approach 
has  been  amply  demonstrated  in  such  applications  as  nonconservative  stabil- 
ity, damping  effects  and  very  general  boundary  conditions  [7,8,9].  In  the 
present  study,  the  solution  formulation  is  limited  to  the  blade  vibration 
considering  only  the  coupling  of  flapping  and  root  torsion.  Although  the 
present  method  can  conceivably  be  extended  for  solutions  of  nonlinear  prob- 
lems, it  is  desirable  first  to  have  a thorough  understanding  to  solutions 
of  linear  problems.  For  this  purpose,  one  can  go  back  two  decades  and  use 
the  equations  consistently  derived  by  Houbolt  and  Brooks  in  1956  [10].  The 
original  set  of  equations  was  derived  for  elastic  distributed  torsion.  It 
can  be  adapted  to  model  root  torsion  if  proper  boundary  conditions  are 
introduced.  This  is  shown  in  Sections  1 and  2.  The  physical  justification 
for  emphasizing  root  torsion  over  the  distributed  torsion  was  due  to  the 
pitch  control  link  at  the  inboard  end  of  the  blade  and  was  used  by  Miller 
and  Ellis  [3]  and  again  by  Ham  [4].  It  should  be  clear  that,  in  the 
present  formulation,  to  include  distributed  torsion  is  simply  a matter  of 
increasing  the  number  of  degrees-of-freedom  of  the  discrete  system.  The 
basis  of  the  solution  formulation  and  the  technique  of  handling  the  coupling 
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terms  are  given  in  Sections  2 and  3.  Finally,  the  numerical  results 
obtained  indicate  that  "flutter  instability"  can  occur  simply  due  to  the 
coupling  effect  considered  without  any  aerodynamic  loads. 


2.  STATEMENT  OF  THE  PROBLEM  - DIFFERENTIAL  EQUATIONS.  As  a first 
step  to  demonstrate  the  application  of  the  unconstrained  variational  - 
finite  element  formulation  to  helicopter  rotors,  the  vibration  of  a rotor 
blade  considering  the  coupling  of  flap  and  root  torsion  modes  of  motion 
is  analyzed  (Figure  1).  For  this  purpose,  the  linear  set  of  equations, 
derived  by  Houbolt  and  Brooks  [10],  including  the  coupling  flap  and  distri- 
buted torsion  is  rewritten  here: 


(EIw'1  - TeA<J>)"  - (Tw')'  - (fllmxe$)»  + m(w  + e<J>)  = lz,  (1) 

- [(GJ  + TkA2)<t>'] ' - TeAw"  + fi2mxew' 

+ ft2m(k2  - k2  - + ee„)d>  + mk2  <p  vineii  = M (2) 

*■  m2  ml  o'*  m Y 

where  w = w(x,t)  and  <J)  = <))(x,t)  denote  the  flapping  deflection  and  distri- 
buted torsion  of  the  rotor  blade  respectively.  A prime  (')  denotes  differ- 
entiation with  respect  to  x,  the  coordinate  along  the  blade  elastic  axis 
and  a dot  (•)>  differentiation  with  respect  to  the  time  t.  Other  symbols 
are  defined  in  the  following  and  are  consistent  with  the  notation  in 
reference  [10]. 

El  = flexural  rigidity  of  the  cross  section. 

GJ  = torsional  rigidity  of  the  cross  section. 

£,  = length  of  the  blade. 

e = chordwise  distance  between  elastic  axis  (E.A.)  and  centre 

of  gravity  (C.G.)  of  a cross-section,  positive  if  C.G.  is 
ahead  of  E.A. 

eA  = chordwise  distance  between  E.A.  and  centre  of  tensile  area 
(C.T.)  of  a cross  section,  positive  if  C.T.  is  ahead  of  E.A. 

eD  = chordwise  distance  at  the  root  between  E.A.  and  the  axis 

about  which  the  blade  is  rotating,  positive  if  E.A.  is  ahead. 


= polar  radius  of  a gyration  of  the  tensile  area  of  a cross- 
section  w.r.t.  E.A.  in  a cross-section. 

= radius  of  gyration  of  the  total  area  of  a cross-section  about 
the  major  neutral  axis  (axis  1-1  in  Figure  2). 

= radius  of  gyration  of  the  total  area  of  a cross-section  about 
an  axis  perpendicular  to  the  major  neutral  axis  and  through 
E.A.  (axis  m2-m2  in  Figure  2). 


k = polar  radius  of  gyration  of  the  total  area  of  cross-section 

about  E.A.  (k2  = k2  + k2  ). 

m ml  m2 

fi  = blade  angular  velocity. 

T = T(x)  = / mftzxdx,  which  is  the  tensile  force  at  location  x. 
x 

Lz  = aerodynamic  lift  per  unit  length  of  the  blade. 

M = aerodynamic  torque  per  unit  length  of  the  blade. 

In  reference  [10],  a more  general  set  of  linear  equations  with  fully 
coupled  flap,  lag  and  distributed  torsion  modes  of  motion  were  derived. 

The  coupling  terms  have  been  shown  to  be  due  to  the  noncoincidence  of  the 
elastic  centre,  centre  of  gravity  and  tension  centre,  the  centrifugal  force 
and  the  built-in  angle  of  twist.  For  a rotor  blade  without  the  built-in 
twist  angle,  the  lag  mode  of  motion  is  uncoupled  from  the  general  set  of 
equations  and  the  remaining  coupled  equations  are  Eqs.  (1)  and  (2) 
considered  here. 

Since  only  the  "free  vibration"  of  the  rotor  blade  with  a coupling 
between  flexural  and  root  torsion  is  considered  in  this  paper,  the  terms 
due  to  aerodynamic  forces  are  set  to  zero  and  the  torsional  displacement 
is  only  a function  of  time  and  not  a function  of  x.  Thus 

M = Lz  = 0,  (})  = $(t)  (3) 

and  Eqs.  (1)  and  (2)  reduce  to  the  following: 

EIw""  - (Tw')'  + mw  - Sl2me<f>  + me$  = 0 

fi2m(k2  - k2  + ee  )<)>  + mk2  <£  - Te  w" 
m2  ml  ° m A 

+ mfi2exw'  + mew  = 0 

In  Eqs.  (4)  and  (5),  it  is  also  assumed' that  the  blade  has  constant 
E,  I and  m throughout  its  length.  To  simplfy  solution  formulations  as 
much  as  possible,  Eqs.  (4)  and  (5)  will  be  transformed  into  dimensionless 
forms  and  appropriate  dimensionless  parameters  will  be  introduced.  This 
process  will  also  facilitate  parametric  studies. 

Let 


w - J . * = t 


where  the  constant  c has  a real  time  dimension  and  will  be  defined  later  in 
Eqs.  (12).  Eqs.  (4)  and  (5)  become 


(4) 

(5) 
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A - A A (T  ♦ A - mrt?e4)  -*■  ~ A|  = 0 
* ai4  £4  ax  ax  c2  at2  c2  at2 


mft2k2<J> 


- -4  (T 

4 21 

3x  3x 

c‘ 

a2$  TeA£ 

32  w 

at2  " £2 

8x2 

ax  c2  at2 


where 


k2  = km2  " kml  + eeo 


Multiplying  Eq.  (7)  by  t one  has 

El 


a“w  _3_  T£2  3w  m£^  1 Aw  m£^  e l 2 

3i4  " ax  ei  ai5  + ei  c2  at2  ei  z c2  c * 

m£4  1 e az4> 

+ r rr  = 0 

ei  c2  z at2 


Multiplying  Eq.  (8)  by  -p-  , one  has 


k2  , , k m a2$  c2T  eA  a2w  . , , e . 3w  , e a2w 

c n z2  $ z2  at2  m£2  £ ax2  c Q £ x ax  £ at2  = 0 


Now  let 


,2  _ m£2 


6 = I ' eA  = X 


r k r Km 

k K r . K " — 


n = cn 


T(i)  = T-(x-I£  - = 1 ft 
EI  2 


ri2(i-x2) 


Eqs.  (10)  and  (11)  then  become 

A*  _ ± (f  4 A 

3x4  3x  3x  3t 


j - 4 (TOO-3*)  * - a’ H * i -3't-  ■ 0 


4 


and 


fi2k2*  + km  7ii-T00  eA^+^ex^.e 


32<J> 


at 


r rZ\  ~ 3 w A o2"C  3w  A - a W 

a2t2 


ax" 


= o 


ax 


(14) 


With  all  quantities  in  Eqs.  (13)  and  (14)  in  dimensionless  forms, 
one  can  omit  the  bars  altogether  and  write: 

w”"  - (Tw')'  + w - n2e<f>  + 4’  = 0 (IS) 

and 

fi2k2<j>  + k2  $ - TeAw"  + $l2exw'  + ew  = 0 (16) 

Furthermore,  it  is  assumed  that 

w(x,t)  = w(x)eXt  (17) 

d>(t)  = <J)eXt 

Thus  the  final  set  of  equations  upon  which  the  present  solution  formula- 
tions are  based,  is  the  following: 

w"'»  - (Tw1)'  + X2w  - fl2 etji  + X2e<f>  = 0 (18) 

and 

n2k2d>  + A2k2  <J)  - TeAw"  + fi2exw'  + X2ew  = 0 (19) 

3.  AN  UNCONSTRAINED  VARIATIONAL  STATEMENT  AND  BOUNDARY  CONDITIONS. 

Some  of  the  unique  features  of  the  unconventional  variational  formulation 
are  that  all  the  boundary  conditions  are  natural  boundary  conditions  and 
that  the  set  of  the  differential  equations,  together  with  all  the  boundary 
conditions,  is  the  direct  consequence  of  a variational  statement  and  vice 
versa.  The  construction  of  such  a variational  statement  is  simply  a process 
of  integration-by-parts  from  a bilinear  functional  of  the  original  differ- 
ential equations  multiplied  by  the  variation  of  the  adjoint  field  variable, 
into  some  other  bilinear  functional  with  lowest  possible  derivations  of 
both  the  original  field  variable  and  the  variations  of  the  adjoint  variable. 
Kith  a proper  choice  of  the  generalized  Lagrange  multipliers,  any  physical 
meaningful  boundary  conditions  consistent  with  the  physical  meaning  of  the 
given  differential  equations  themselves  can  be  resulted  from  the  variational 
statement  as  natural  boundary  conditions.  This  general  process  was  treated 
elsewhere  [7]  and  will  not  be  reported  here.  Presently,  an  unconstrained 
variational  statement  will  be  given  which  leads  to  the  original  differential 
equations  and  a set  of  a very  general  boundary  conditions. 
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Let  us  consider  the  variational  statement 

61  = 0 (20a) 

with 

I = Ij  + I2  (20b) 

1 1 = J [W'w*"  + Tw'w*'  + A2ww*  + e(A2  - ft2)4>w*]dx 
0 

+ h1w(0)w*(0)  + k2w>  (0)w*‘  (0)  (20c) 

and 

1 2 ~ I t^2k-2  + ^2k^)<Jxl>*  + eAT'w'(J)*  + e(Q2xw'<|)*  + X2w$*)]dx  + k, <}>()>* 

0 (20d) 

where  a star(*)  denotes  the  adjoint  variables  and  the  variations  are  totally 
unconstrained.  The  Lagrange  multipliers  kj,  k2  and  kj  are  the  spring  con- 
stants, in  dimensionless  form,  for  deflection,  bending  and  torsion  at  the 
hub  (x  = 0)  respectively.  To  show  that  the  unconstrained  variational  state- 
ment of  Eqs.  (20)  leads  to  the  original  differential  equations  and  the 
necessary  boundary  conditions,  one  considers 


(61)  , = 0 

w,<p 


(21a) 


where  (6I)w>(j)  means  taking  the  variation  of  I with  ui,<J>  not  varied.  Thus 

(61)  = (610  + (6I2)  =0  • (21b) 

w,<P  w, q>  ^ w,<j> 


with 


(6lj)w  ^ = /1[w"5w*"  + Tw'6w*'  + A2w6w*  + e(A2  - n2)<}>6w*]dx 


and 


+ k1w(0)6w* (0)  + k2w,'(0)6w*' (0) 

(6i2)w,$  = + l2^*6**  + eAT,w'6^* 

+ e(n2xw'  + A2w)6(J)  ]dx  + kj<})6<}>* 


(21c) 


(21  d) 
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Performing  integration-by-parts , one  arrives  at 

(<5li)  a = /1[w""  - (Tw')’  + A2w  + e(A2  - n2)<J>]6w*dx 
w»  0 

+ w"(l)6w*'  (1)  - [w"tl)  - T(l)w'  (l)]6w*(l) 

- [w"(0)  - k2w' (0)]6w*‘ (0) 

♦ [«’"  (0)  - T(0)w'(0)  + k1w(0)]6w*(0)  (22a) 


and 


(6i2)  = / [(fi2k2  + A2k2)<J>  - e.Tw"  + e(n2xw*  + A2w)  ] 6<J>*dx 

W»<P  o m a 

+ [eAT(l)w'(l)  - eAT(0)w'(0)  + k3<J>] 6<J>*  (22b) 

Thus  Eqs.  (21)  is  the  necessary  and  sufficient  condition  for  the  following 
differential  equations  and  boundary  conditions.- 


D. E. : 


w""  - (Tw')'  + A w ♦ e (A2  - ) 4>  = 0 


(fi2k2  + A2k2)<J>  - eATw"  + e(fi2xw’  + X2w)  = 0 


B.C. : 


and 


w"(l)  = 0,  W"  (1)  = 0 
w"(0)  - k2w'  (0)  = 0 


o2 

w"'  (0)  - y-  w' (0)  + klW(0)  = 0 


e. 

- n2w  (o)  + k3<j)  = o 


(23a) 

(23b) 

(24a, 24b) 
(24c) 

(24d) 

(24e) 


Note  that  in  Eqs.  (24)  the  fact  that 

T(1 ) = 0 

and 


T (0)  = I n2 


(25) 


has  been  used. 


a 


■ — ~ — — - — m i i — 


It  is  observed  that  the  boundary  conditions  of  Eqs.  (24c)  and  (24d) 
represent  very  general  support  conditions  at  the  hub.  The  special  case 
kj  = «>  and  k2  = 0 corresponds  to  a hinged  blade,  while  the  case  kj  = ®, 
k2  = 00  corresponds  to  a hingeless  blade.  These  two  special  cases  are  the 
only  ones  considered  in  the  literature.  The  boundary  condition  (24e) 
indicates  a coupling  between  the  flexural  motion  and  root  torsion  at 
the  hub  due  to  the  centrifugal  forces  and  the  noncoincidence  of  the  elastic 
axis  and  the  tension  axis. 

4.  MATRIX  EQUATIONS  FROM  FINITE  ELEMENT  DISCRETIZATION.  In  this 
section,  we  shall  briefly  describe  the  formulation  of  the  matrix  equation 
of  the  approximate  solution  from  the  variational  statement  given  in  the 
previous  section.  From  Eqs.  (21),  one  can  write 


C«)w  * = 0 ; 

w,  <J> 


(61)  = /1(w"6w*"  + Tw'6w*'  - efi2<fi6w*)dx 

0 

+ klW(0)6w*(0)  + k7w’  (0)6w*’  (0)  + X2  /*  (w6w*  + e<j>6w*)dx 

0 


+ /*  [fi2k2<f>6<f>*  + eAT'w'6(J)*  + efi2xw' 6<J)*]dx  + k-j(J>6<J>* 

0 

+ X2  (k2<}>6<}>*  + ew6(J>  )dx  (26) 

0 m 

It  should  be  noted  that  Eq.  (26)  is  a quite  general  equation.  Various 
types  of  approximate  solutions  can  be  obtained  depending  on  the  choice  of 
the  coordinate  functions.  The  motivation  of  using  the  finite  element 
discretizations,  which  corresponds  to  the  choice  of  a set  of  piecewise 
analytic  functions,  is  twofold:  for  the  ease  of  extending  this  formulation 
to  problems  of  irregular  geometry  and  for  the  adaptations  to  general  finite 
element  computer  systems. 


In  the  present  formulation,  however,  only  blades  of  uniform  cross- 
sections  will  be  considered  and  the  elements  are  assumed  to  be  of.  the  same 
length.  Thus,  one  introduces  a local  (element)  coordinate  £ which  relates 
to  the  global  (entire  blade)  coordinate  x such  that 


K = €(i)  = L(x  - i=i)  (27) 

L 

where  L is  the  number  of  elements  and  i denotes  the  i-th  element.  Using 
the  notation 

w(i)(0  = w[x(C(l))],  etc.  (28) 

for  simplicity  and  noting  that 

d£  = Ldx 


w'(x)  = Lw(l),(£),  etc. 


(29) 
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one  then  has,  from  Eq.  (26): 


(61)  =0 
w,4> 

L i 

= I {/  L3w(i)"(C)6w*(i)"(£)d£ 

i=l  0 

t /XL  T(l)(C)w(l),(5)6w*'(l)(adC  - ^4^  /16w*(i5d5) 
0 . 0 

+ (0)6w*t^  (0)  + k2L2wtl)'  (0)6w*(l^'  (0) 

+ X2  £ {I  /1w(i)(06w*(i)(5)d^  + /16w*(i)(C)d?} 

i=l  L 0 L 0 

-l  {n2k2<J.6<{.*  Y"  /1dCCl^  + eA  L f\™'  (C)w(l)'(C)dC6$* 


+ "J1^  + (i  - l)]wCO,(€)de  64.*}  + k34>64>J 


♦ x2  { —■  4)64>"  /Vi5  * l f r« 

L 0 i=lL  0 


rl  (i) 


K)dU«* 


where 


T^U)  = - b[^z  + 2 (i  - 1)C  + (i  - D2  * k2] 


At  this  point,  it  is  appropriate  to  introduce  the  shape  function  and 
generalized  coordinates.  Let 

w(l)(£)  = aT(C)W(i) 


w*(i)K)  = aT(C)W*fi5 

where  the  superscript  T denotes  the  transpose  of  a matrix, 
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and 


»WT  = { 

w2^ 

»3(i’ 

W4(i)  } 

w*(i)T  _ { w*(i) 

»;(i) 

»;(i) 

) 

(34) 

352  + 2?3,  ? - 

2^2  ♦ 5 

\ 3£2 

- 25’,  - e * Cs) 

(35) 

Using  Eqs . (31),  (32)  and  in  terms  of  W^,  Eq.  (30)  becomes 


(«)„  A = o 

w,<f» 


= l 6W*(1)T{L3/1a"a"Td£  - bL[/1C2a* a,Td£  + 2 (i  - l)f\ a'a,TdC 
i-1  0 0 0 


{L2  - (i  - l)2}  /1a,a,Td5]W^1^} 
0 

r *(1)T  rl 

+ l 6W  1 ' (-  2ebL)  / ad£  $ 

i=l  0" 


+ 6W*  t1)1  kx  a(0)aT(0)  W*-1^  + k2L2a’  (0)a,T(0) 


+ X2  { £ 6W  ^ j-  f aaTd£  + j-  J ad£  <J>) 
i=l  ~ 0 0 ~ 

+ 6<f>  f22kz4>  + 6<J>  kj  ij) 

+ 6<f>*  l 2 (e  - e )bL  /V'T(£)d£  + (i  - l)/^'1^  W 


(i) 


i=l 


A2  6<t>*[k2<j)  + F f f1aTdE>  W(l)] 
m i=l  L 0~  ~ 


(36) 
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I 


(«I)  = 0 

w > <p 


= l 6W*UJ  {L3C  - bL[E  + 2 (i  - 1)D  - (L2  - (i  - l)2)  B] } W(l) 
i=l  ~ ~ ~ 

+ l SVi*^T  (-  2 ebL)p  <}>  + 6W*(i)T  [k.F  + k,L  GlW^1^ 
i=l  ~ - 1 ~ - 

+ \2[  l 6W*^)T  [ i A W(l)  ♦ f p <fe] 
i=l  L ~ ~ L - 

+ 6.J.*  l 2 (e  - e.)bL  (sT  + (i  - l)rT)  W(l) 
i=l  ~ 

+ 6<J>  (ft2k2  + k,)<J>  + A2  64>*  ( l f pT  + k2d>) 

i=l  L ~ m 


where  the  element  matrices  are  defined  as  follows, 
A = /^aa^rir.  r = r 


A = fxa.ald%,  B = f1a'a'  *d£,  C = /1a"a"Td£ 

o—  o - o~  - 

D = /1Ca'a,Td5,  E = /Va'a'Td£ 


F = a(0)  aT(0) , G = a' (0)a'T(0) 


,1  -1 
p = / a dC,  r * / a'  d£, 

“*  A ^ A ^ 


s = / 5 a'  d£ 
0 ~ 


0~  0 ~ 0 ~ 

The  numerical  values  of  these  matrices  are  given  in  the  Appendix. 
In  terms  of  the  global  coordinates  defined  in  the  following 


-“!;(•  Mil 


II 

rS 

(i) 

H 

W 

CM 

»3(1) 

»4(i) 

W3(a 

K (2) 

"4 

W (L) 

• -w3 

W 

W4 

(i) 

w *0) 
W2 

*;(1) 

»;tl) 

„*(2) 

*3 

w*(2) 

"4 

• -"3 

W*(L) 

*'4 

I 


r 


One  has  from  Eq.  (37): 


(61>w,*  = 0 


= 6W  1 (K  + X2M)W 


{6W*T64>*} 


K11  + * M11  K12  + * M12 


K21  + **M21  K22  + * M22 


W 


(41) 


The  global  matrices  K^j,  Mij,  i,j  =1,2  are  assembled  from  the  element 
matrices  defined  in  Eqs.  (38)  in  the  following  manner. 


6W*TKnW  = l 6W*UJ1  (L3(C  + bB)  + bL(2D  - B - E) 


(i)T 


i=l 


+ 2ibL(B  - D)  - i2bLB)W(l-)  + 6W*(1)T{k1F  + k2L2G)K(1^ 


■*Ci)T 


6W*TK124>  = l 6W*^J'(-  2ebL  p)<f> 


(42a) 

(42b) 


i=l 


J , ^T1UI(i) 


6<f>*K?1  W = 6$*  l 2 (e  - e.)bL[sT  - r + ir']w' 
i=l 


6<f>*K22<t>  = 6<j>*(fi2k2  + k3)<f> 
*M  .W  = l 6W*(l)  y:  A W(l) 


6W  ■•■j  ji. 


i=l 


(42c) 

(42d) 

(42e) 


6W*M.  74»  = l f p <{> 

~ i=l  L ~ 


(42f) 


6<f*M 91W  = 6$*  l ± pT  W(l) 
i=l  b ~ ~ 


21~ 


6<f>  M 22<p  = 6<p  k2  4> 


(42g) 

(42h) 
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r-~ — 


r 


Since  6W  ^ = (6W  ^ 6<f>  } in  Eq.  (41)  is  unconstrained,  Eq. 
directly  to 

(K  ♦ X2M)W  = 0 

which  is  the  final  matrix  eigenvalue  equation  to  be  solved. 

S.  RESULTS  AND  DISCUSSION.  Prior  to  the  presentation  of 
strative  numerical  results,  it  will  be  worthwhile  to  maVe  some 
on  the  nondimensionalized  differential  equations  (15)  and  (16) 
boundary  conditions  (24). 

a.  For  e = 0,  the  flapping  motion  w(x,t)  and  the  root  torsion  <|>(t) 
are  essentially  uncoupled.  The  eigenvalue  solutions  of  w reduce  to  that  of 
a rotating  beam  and  agree  well  with  the  available  data  of  Boyce,  DiPrima 
and  Handelman  [11].  The  torsional  vibration  frequency  has  only  one  eigen- 
value solution  and,  if  ey^  = 0 also,  it  varies  linearly  with  the  rotor  speed 
ft  as  expected. 

b.  For  e ^ 0 and  e^  = 0,  the  motions  are  generally  coupled.  It  is 
observed  that  if  e and  <J>  both  change  sign,  the  equation  remains  unaltered. 
Thus,  as  far  as  eigenvalues  are  concerned,  they  depend  only  on  the  absolute 
values  of  e.  The  solution  of  <f>  for  a negative  e,  however,  is  180°  out  of 
phase  compared  with  the  one  for  a positive  e of  the  same  magnitude. 

c.  For  e / 0 and  e^  i 0 the  motions  are  generally  coupled.  It  is 
observed  that  if  e,  e.  and  <{>  all  three  change  signs,  the  governing  equa- 
tions remain  unaltered. 

Some  demonstrative  calculations  will  now  be  given*.  The  eigenvalue 
X is  generally  a complex  number. 

X = XR  + iXj  (44) 

From  equations  (17)  and  (44),  it  is  clear  that  the  system  is  unstable  of 
divergence  when  Xr  is  nonzero  positive  and  Xj  = 0.  When  Xr>  = 0,  on  the 
other  hand,  Xj  (LAMBDA)**  represents  the  nondimensional  frequency  of 
vibration  and  it  can  then  be  plotted  against  the  nondimensional  blade 
rotating  speed  ft  (OMEGA),  as  shown  in  Figures  3 through  7.  When  X is 
complex,  one  of  the  square  roots  of  X2  must  have  nonzero  positive  real 
part.  The  system  is  then  unstable  since  X appears  in  the  equations  only 
as  X2.  The  value  of  X generally  become  complex  as  the  two  branches  of 
the  frequency  curve  coalesce.  The  "critical"  speeds  can  be  located  in 
these  figures  by  noting  these  points  of  coalescence. 


*An  extensive  parametric  study  of  a rotor  blade  instability  in  vacuum  due 
to  the  coupled  flap- (root)torsion  motion  will  be  presented  in  a 
forthcoming  paper. 

**The  symbols  enclosed  in  parentheses  are  those  used  in  Figures  3 
through  7. 


(41)  leads 


(43) 


the  demon - 
observations 
and  the 
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For  the  following  sample  calculations,  the  parameters  used  are  typical 
for  a medium  size  rotor  blade  of  20  ft.  in  length,  for  example.  It  is 
further  assumed  that  (the  nondimensionalized  quantities) 

e (E)  = 0.003,  e0(E0)  = 0 

kml(KMl)  = 0.0025,  km2  (KW2)  = 0.0100 

The  boundary  conditions  at  the  hub  are  those  of  a hingeless  blade.  The 
values  of  kj(K-l)  and  k2(K-2)  are  set  to  108  as  approximations  to  infinity. 
The  torsional  spring  constant  kj(K-3)  has  been  set  to  zero. 

In  Figure  3,  e^(E^)  is  taken  to  be  0.003.  Thus 

y (GAMMA)  = ^ ~ = 1.00  . 

Here  the  lowest  branch  of  Aj  is  essentially  for  the  torsional  motion.  The 
second  and  third  lowest  branches  are  essentially  the  first  and  the  second 
for  flapping  motion.  For  the  range  of  rotor  speed  shown*,  0 <_  n <_  25,  the 
coupling  is  not  sufficient  to  have  instability.  In  the  subsequent  figures 
as  e^  = 0.0015,  0,  -0.0015  and  -0.003  (and  Y = 0.5,  0,  -0.5  and  -1.0), 
the  effect  of  coupling  becomes  more  and  more  severe.  In  Figure  4,  the 
two  lowest  branches  of  eigenvalues  appear  to  draw  closer  compared  with 
those  in  Figure  3.  They  actually  coalesce  in  Figure  5 at  a critical 
speed  about  fl  = 16.5.  As  e^  continue  to  decrease  (increase)  in  algebraic 
sense  while  holding  e a positive  (negative)  constant  the  critical  speed 
of  flutter  instability  tends  to  decrease  and  thus  the  structure  becomes 
more  critical.  This  is  observed  in  Figures  6 and  7. 

To  obtain  some  idea  on  the  effect  due  to  the  torsional  spring  k , the 
calculations  are  repeated  for  the  same  set  of  data  as  in  Figure  5 except 
that  the  value  of  k^  is  changed  from  zero  to  0.01.  The  new  frequency  curves 
are  shown  in  Figure  8.  It  appears  in  this  figure  that  the  essentially 
torsional  branch  has  moved  upward  while  the  other  branches  remain  largely 
unaltered.  The  new  critical  speed  prior  to  flutter  instability,  however, 
has  changed  from  16.46  (Figure  5)  to  17.08. 

A rather  quick  review  of  the  existing  literature  suggests  that  the 
flutter  instability  due  to  the  coupling  flapping  and  root  torsion  motions 
may  have  not  been  reported  previously.  The  practical  significance  of  such 
instability  within  the  range  of  real  design  parameters  is  obvious.  An 
extensive  parametric  study  on  this  instability  is  being  carried  out.  The 
results  will  be  reported  shortly. 
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*At  A (OMEGA)  = 25,  which  is  nondimensional , the  real  rotor  speed  is  roughly 
1000  rpm  for  a typical  value  of  c = /(pAJlVEI)  = 0.0040min. 
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APPENDIX 


The  numerical  value  of  some  of  the  matrices  used  in  Section  4 are 
given  here. 


A = / aa  a£ 
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B = / a’a’  d£  = 
0 ~ ~ 
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a"a"  d£  = 


12  6 -12  6 

6 4 - 6 2 


m 


T.A.  = Axis  of  Torsion 
E.A.  = Elastic  Axis 
C.G.  = Centre  of  Gravity 
C.T.  = Centroid  of  Tensile  Area  of  a Cross  Section 


Figure  2.  Parameters  related  to  off-sets  of  various  axes. 
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Figure  3.  Vibration  frequency  vs  speed  of  rotation  (eA*0.0030,  e=0.0030). 
No  instability  in  the  range  of  speed  shown. 
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Figure  4.  Vibration  frequency  vs  speed  of  rotation  (eA=0.0015 
No  instability  in  the  range  of  speed  shown. 
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Figure  5.  Vibration  frequency  vs  speed  of  rotation  (eA=0,  e=0.0030). 
Flutter  instability  at  (M6.46. 
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Figure  6.  Vibration  frequency  vs  speed  of  rotation  (eA=-0.0015,  e=0.0030). 
Flutter  instability  at  0=5.80. 
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Figure  7.  Vibration  frequency  vs  speed  of  rotation  (eA«-0.0030,  e*0.0030). 
Flutter  instability  at  ft»4.59. 
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Figure  8.  Vibration  frequency  vs  speed  of  rotation.  (Same  data  as  in 


Figure  5,  except  IC3  has  been  changed  from  zero  to  0.01). 
Flutter  instability  at  fi=17.08. 
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